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ABSTRACT 


The  boundary  layer  thickness  around  a  bLunt -nosed  body 
increases  as  the  Reynolds  number  decreases.  When  this  thickness 
reaches  a  magnitude  of  the  order  of  the  inviscid  shock  detachment 
distance,  the  heat  fluxes  and  the  components  of  stresses  after  the  shock 
begin  to  play  important  roles  in  the  fLow  field  behind  the  shock.  These 
effects  necessitate  the  modification  of  the  usual  Rankine-Hugoniot  shock 
conditions  in  the  study  of  the  energy  transfer  between  the  fluid  and  the 
body  surface. 

In  the  present  anaLysis  the  leading  order  shock  conditions  are 
obtained  from  the  systematic  order  of  magnitude  analysis.  A  complete 
matching  scheme  for  the  components  of  heat  flux  and  stresses  between 
the  shock  and  the  boundary  layer  is  proposed  for  the  study  of  the  flow 
field  as  well  as  the  surface  heat  transfer  to  the  body.  No  postulation, 
except  the  axial  symmetry,  has  been  made  concerning  the  shape  of  the 
shock.  The  shape  of  the  matching  surface  can  be  obtained  uniquely 
from  the  present  analysis.  A  method  of  solution  has  been  introduced  so 
that  the  laborious  job  of  numerical  integration  of  the  flow  equations  with 
complicated  boundary  conditions  can  be  bypassed.  The  governing  equa¬ 
tions  are  transformed  and  reduced  to  algebraic  forms  with  the  knowledge 
of  the  usual  similarity  solutions  of  the  boundary  layer  equations. 

Calculation  has  been  performed  for  the  specific  case  of  M0O=5.7 
and  T0  =  the  free  stream  stagnation  temperature  =  2100°R.  The  heat 
transfer  result  at  the  stagnation  point  is  in  good  agreement  with  the 
experimental  value  and  shows  a  slight  improvement  to  Cheng’s  Thin 
Shock  Layer  Theory  in  the  range  of  application  of  the  present  schemes. 
The  density  and  the  stagnation  enthalpy  at  the  matching  point  are  obtained 
as  functions  of  the  Reynolds  number.  They  both  show  the  appropriate 
asymptotic  values  of  a  discontinuous  shock  as  the  Reynolds  number 
increases.  The  matching  distance  and  the  shock  eccentricity  are  also 
obtained  as  functions  of  the  Reynolds  number. 
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SECTION  I 


INTRODUCTION 

At  high  Reynolds  number,  the  problem  of  hypersonic  cold  blunt 
body  heat  transfer  is  considered  to  be  solved  within  the  framework  of 
the  boundary  layer  analysis  with  negligible  boundary  layer  thickness'  *  ' . 
The  thin  boundary  layer  near  the  wall  remains  to  be  a  good  approximation 
to  the  flow  field  around  the  body.  Consistent  with  the  cold  wall  and  the 
hypersonic  assumptions,  similarity  solution  of  the  boundary  layer  equa¬ 
tions  can  be  obtained  with  the  usual  classical  boundary  conditions.  As 
the  free  stream  density  decreases,  this  model  no  longer  provides  an 
adequate  picture  of  the  flow  field  around  the  body,  and  it  is  necessary 
to  study  the  errors  imbedded  in  the  scheme. 

The  effect  of  the  vorticity  interaction  and  its  importance  to  the 

flow  characteristics  around  the  body  has  been  first  pointed  out  by  Ferri 
(3) 

and  Libby'  .  The  shock  generated  vorticity  in  the  inviscid  flow  behind 
the  shock  is  a  function  of  the  detachment  distance,  and  increases  with 
an  increase  of  free  stream  velocity.  On  the  other  hand,  the  average 
shearing  stress  in  the  boundary  layer  decreases  as  a  result  of  increasing 
in  the  boundary  layer  thickness.  When  the  shearing  stresses  in  the 
boundary  layer  and  the  inviscid  shear  layer  outside  become  comparable 
to  each  other,  the  interference  effect  which  exists  between  the  two  layers 
becomes  important  to  the  flow  field  behind  the  shock.  Investigation  has 
been  made  on  this  effect,  and  theory  has  been  developed  in  reference  4 
on  the  heat  transfer  of  a  blunt  hypersonic  geometric  configuration  as  a 
function  of  parameters  characterizing  the  vorticity  and  the  boundary 
layer  thickness  after  the  shock.  The  theory  is  based  on  a  postulation 
that  for  a  cold  blunt  body  the  boundary  layer  displacement  effect  is 
negligibly  small  as  compared  with  the  effect  of  the  vorticity  interaction. 
With  a  boundary  layer-inviscid  shear  layer -discontinuous  shock  model. 
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the  boundary  layer  solution  has  been  obtained  by  matching  both  the 
velocity  and  the  shear  with  the  inviscid  layer  at  a  point  determined  by 
matching  the  mass  flow.  The  experimental  data  on  heat  transfer 
value^  show  that  the  theory  is  verified  down  to  the  range  of  Reynolds 
number  based  on  the  quantities  after  the  shock  of  the  order  of  500  or 
more.  At  this  range  of  Reynolds  number,  an  estimation  shows  that  the 
boundary  layer  thickness  is  almost  of  the  order  of  the  inviscid 
detachment  distance.  The  theory  breaks  dpwn  beyond  this  point  where 
the  stress  components  and  the  components  of  the  heat  flux  after  the 
shock  begin  to  play  an  important  role  in  the  flow  field  around  the  body^. 

Van  Dyke^  has  made  the  general  discussion  on  the  second  order 
effects,  such  as  the  curvature  effect,  the  velocity  slip  and  temperature 
jump  effects  at  the  wall,  the  vorticity  effect  and  the  displacement  effect 
due  to  the  first  order  boundary  layer.  Using  the  usual  inner-outer  two 
layer  expansion  scheme,  the  fir3t  order  and  the  second  order  equations 
are  presented  with  boundary  conditions  obtained  by  using  the  standard 
asymptotic  matching  technique.  The  second  order  effects  are  expressed 
in  terms  of  the  first  and  the  second  order  solutions.  In  the  formulation 
of  the  problem,  the  outer  layer  as  defined  in  reference  7  has  been 
treated  without  considering  the  existence  of  the  shock  wave.  This  implies 
that  the  range  of  Reynolds  number  which  the  theory  is  to  apply  has  to  be 
high  to  ensure  the  convergence  of  such  an  expansion  scheme.  The 
boundary  layer  thickness  is  much  smaller  than  the  inviscid  detachment 
distance.  At  this  range  of  Reynolds  number,  the  effect  of  the  dis¬ 
placement  thickness  may  be  of  the  same  degree  of  importance  as  the 
effect  of  the  vorticity  interaction.  It  is  known  that  to  obtain  the  first 
order  boundary  layer  solution  requires  the  knowledge  of  the  first  order 
inviscid  solution,  whiLe  the  second  order  boundary  layer  solution 
requires  the  knowledge  of  the  next  order  solution  in  the  outer  layer. 

This  theoretically  speaking,  can  be  obtained  only  by  solving  the  boundary 
value  problem  between  the  perturbed  body  and  the  perturbed  shock,  if 
the  displacement  effect  is  to  be  considered  consistently.  In  the  appli¬ 
cation  to  the  blunt  body  problem  in  reference  7,  numerical  results  have 
been  presented  for  all  the  second  order  effects  except,  the  displacement 
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effect,  due  to  the  extreme  complication  involved.  No  new  scheme  has  been 
proposed  concerning  this  matter.  The  effects  of  velocity  slip  and  temper¬ 
ature  jump  reported  in  reference  7  have  been  questioned  by  Cheng  in  a 
recent  review^ ®),  since  these  effects  cannot  be  important  insomuch  as  the 
problem  of  cold  blunt  body  in  hypersonic  flow  is  considered. 

Probstein^  has  proposed  the  classification  of  the  "flow  regimes" 
according  to  the  degree  of  rarefiedness  of  the  surrounding  flow  field  of  a 
blunt  body  reentering  the  atmosphere  with  hypersonic  speed.  In  a  com¬ 
prehensive  review  concerning  the  subject these  flow  regimes  are 
discussed  in  a  manner  that  shows  the  general  picture  of  the  development 
of  the  shock  wave  and  of  the  flow  field  as  the  free  stream  density 
decreases.  No  rigorous  treatment  has  been  presented  concerning  the 
ranges  of  applications  of  these  regimes,  therefore  they  cannot  be  regarded 
as  a  conclusive  statement  on  the  subject.  In  fact,  the  existence  of  the 
"Viscous  Layer  Regime"  defined  by  Probstein  has  been  questioned  by 
several  authors Different  schemes  have  been  proposed  for  the  solutions 
in  the  "Vorticity  Interaction  Regime,"  the  "Viscous  Layer  Regime,"  and 
the  "Incipient  Merged  Layer  Regime"  as  defined  by  Probstein.  In 

reference  10,  Probstein  and  Kemp  have  adopted  the  constant  density  model 
in  treating  the  flow  in  the  "Viscous  Layer  Regime."  Such  a  model,  as 
pointed  out  by  Cheng^^,  describes  more  appropriately  the  physical  problem 
of  an  insulated  body.  Therefore,  one  cannot  expect  an  accurate  prediction 
of  the  heat  transfer  value  from  the  theory.  In  a  more  recent  paper.  Ho  and 
and  Probstein^  have  treated  the  same  flow  regime  with  an  improved 
model.  The  flow  equations  behind  the  shock  have  been  obtained  by  simpli¬ 
fying  the  Navier-Stokes  equations  based  on  the  thin  shock  layer  approach. 

The  regular  Rankine-Hugoniot  shock  conditions  and  the  conditions  at  the 
wall  serve  as  the  boundary  conditions  of  the  problem.  The  heat  transfer 
result  from  the  theory  shows  the  general  trend  of  increase  in  heat  transfer 
value  as  the  Reynolds  number  decreases  but  falls  of  the  order  of  ten  per¬ 
cent  based  on  the  boundary  Layer  value  below  the  experimental  data. 
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The  inclusion  of  the  effects  of  the  viscous  stresses  and  the  heat 
conduction  behind  the  shock  has  been  studied  by  Cheng The  flow 
equations  behind  the  shock  have  been  obtained  by  simplifying  the 
Navier-Stokes  equations  consistent  with  the  thin  shock  layer  approxi¬ 
mation.  Two  different  schemes  have  been  proposed  for  the  solution  of 

.  pmUma  T 

the  problem  depending  on  whether  the  parameter  K  -  e (  is 

o 

large  or  small  as  described  in  reference  12.  In  "Regime  I"  or  the 
"Vorticity-Interaction"  regime  as  defined  by  Cheng,  the  solution  has 
been  obtained  by  integrating  the  flow  equations  satisfying  the  six 
boundary  conditions  one  of  which  determines  the  upper  integration 


limit  ng  (boundary  layer  coordinate).  The  inviscid  theory  has  been 
adopted  in  prescribing  the  boundary  conditions  at  n=  n  g»  and  the 
effects  of  the  viscous  stresses  and  the  heat  conduction  behind  the  shock 


have  not  been  included  here  in  "Regime  I."  The  theory  does  not 
predict  the  actual  detachment  distance  m  a  unique  manner.  The  cor¬ 
rection  term  to  the  inviscid  detachment  distance  presented  in  refer¬ 
ence  12  indicates  the  displacement  effects  of  the  flow  field  behind  the 
shock,  which  can  only  be  obtained  after  performing  the  integration. 
Therefore,  it  is  clear  that  the  scheme  is  reasonably  justified  only  if 
the  displacement  effects  are  smalL,  and  we  find  that  this  is  in  agree¬ 
ment  with  the  postulation  made  in  reference  4.  In  "Regime  II" 
defined  by  Cheng,  the  effects  of  the  viscous  stresses  and  the  heat 
conduction  behind  the  shock  have  been  included  in  the  boundary  con¬ 
ditions.  A  perturbation  type  of  approach  with  parameter  charac¬ 
terizing  the  tangential  pressure  gradient  of  the  flow  has  been  adopted 
for  the  solutions  of  the  problem.  Again,  the  detachment  distance  does 
not  come  into  the  analysis.  Therefore,  the  fact  that  the  detachment 
distance  and  the  flow  behind  the  shock  are  coupled  has  been  ignored. 
Since  essentially  the  scheme  still  adopts  the  concentric  shock  model, 
the  effect  of  the  shock  curvature  cannot  be  observed  from  the  analysis. 
(Notice  that  the  shock  curvature  and  the  detachment  distance  are 
different  parameters.) 
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In  the  present  study,  a  compLete  matching  scheme  for  the  com¬ 
ponents  of  the  stresses  and  the  heat  flux  is  introduced  for  the  study  of 
the  flow  field  behind  the  shock.  The  shock  surface  or  the  matching 
surface  can  be  determined  uniquely  from  the  analysis.  The  detail  of 
the  scheme  is  explained  in  a  later  part  of  this  report. 
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SECTION  n 


DISCUSSION  OF  THE  PROBLEM 

The  problem  here  is  to  investigate  the  fLow  Held  about  a  cold 
blunt  nosed  body  in  hypersonic  low  density  stream.  We  make,  as 
usual,  the  following  justifiable  postulations:  (1)  the  disturbed  region 
of  the  flow  field  is  confined  to  a  distance  much  smaller  than  the  nose 
radius  of  the  body;  (2)  the  flow  density  around  the  body  is  much  larger 
than  the  free  stream  density,  and  the  continuum  formulation  is  still 
adopted  in  the  present  study.  Due  to  the  very  large  stagnation  temper¬ 
ature  to  wall  temperature  ratio,  the  temperature  gradient  in  the 
direction  normal  to  the  wall  is  much  larger  than  that  of  the  circum¬ 
ferential  direction.  The  temperature  jump  at  the  wall  cannot  be 
important  insomuch  as  the  problem  of  a  cold  blunt  body  in  hypersonic 
flow  is  considered.  Therefore  the  usual  thermal  boundary  layer 
remains  a  good  approximation  to  the  temperature  distribution,  provided 
appropriate  boundary  condition  can  be  prescribed  at  the  edge  of  the 
temperature  layer.  The  momentum  equation  of  the  boundary  layer 
equations  is  subject  to  the  error  that  the  pressure  change  acrosB  the 
boundary  layer  has  been  ignored.  The  use  of  boundary  layer  equa¬ 
tions  is  justified  by  the  fact  that  either  the  normal  pressure  gradient 
itself  is  small,  or  the  boundary  layer  is  thin.  Since  we  are  treating 
the  flow  field  in  the  shock  layer  which  has  a  dimension  much  smaller 
than  the  nose  radius,  it  is  reasonable  to  assume  that  a  large  portion 
of  the  pressure  jump  is  absorbed  in  the  shock  transition  zone.  This 
iB  confirmed  by  the  recent  work  by  Cheng'  where  the  pressure 
remains  constant  along  the  axis  of  symmetry  has  been  postulated  in  the 
analysis  of  stagnation  region.  Therefore,  with  this  discussion  in 
mind,  we  may  seek  the  solution  of  the  boundary  layer  equation,  with 
the  shock  conditions  serving  appropriately  as  the  boundary  conditions 
of  the  problem.  The  systematic  analysis  of  obtaining  the  shock  con¬ 
ditions  is  discussed  in  the  details  of  the  present  work. 


SECTION  m 


SHOCK  CONSERVATION  EQUATIONS 

We  assume  that  the  Navier-Stokes  equations  hold  for  the  range 
of  investigation.  The  continuity  equation,  the  equations  of  motion  and 
the  energy  equation  can  be  expressed  in  any  curvilinear  coordinate 
system  suitable  for  the  analysis.  For  an  axially  symmetric  blunt-nosed 
body  at  zero  angle  of  attack,  the  shock  surface  close  to  the  axis  of 
symmetry  can  be  approximated  by  a  sphere  with  the  center  on  the  axis. 

In  this  axisymmetric  spherical  polar  coordinate  system,  the  system  of 
steady  Navier-Stokes  equations,  neglecting  the  effects  of  the  body  force, 
can  be  expressed  as  follows^ 

The  continuity  equation  is 

(pv)  +  2  (pu)  +  cot  0=0  (III-- l  a) 


The  radial  direction  of  momentum  equation  is 
2 


otv  It  +  7  If  -  T-]=  w > +  7  le '"re1  +  7  [2nr.-"„n-nee+nrec‘>t 

(m-ib) 


The  circumferential  direction  of  momentum  equation  is 

p[vl7  +  7  +  (ITr0)+r  (IW  +  7  [(IIee"nnn)c0t  6+3IW 

(m-ic) 

The  energy  equation  is 

r  9h  ,  u  9hn  9p  ,  u  9p,lr9  9T»  9  ,k  9T%  »9T  ,  cot  0  9T.-i  ,  - 

p£v-57  +  7  -F0]=v^f  +  7  ^0  +7^{kr^7)  +  -50  (7  +fc%7  +  ”77"  lFF)1  +  $ 

(m-id) 
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where 


e 


r8 


8 

r 


— 

r  '51' 


and  v  and  u  are  the  veLocity  components  along  the  radial  and  circum¬ 
ferential  directions  in  the  meridian  plane;  and  p,  p,  T,  h,  (i,  \i',  k  are 
the  density,  the  pressure,  the  temperature,  the  specific  enthalpy,  the 
shear  and  the  bult  coefficients  of  viscosity  and  the  coefficient  of  heat 
conduction,  respectively. 

System  of  Eqs.  (HI-l)  is  approximated  by  the  technique  of  the 
order  of  magnitude  analysis  developed  in  reference  14.  Introduce  for 
the  dependent  variables  v,  u,  p,  etc.  the  following  expansion  scheme  in 
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the  shock: 


v(s,0,v)  =  +  vv^  + 


o(v) 


(m-2) 


u(s,0,v)  -  u<°>  +  +  o(v)  etc., 

where  s  is  the  stretched  radial  coordinate  across  the  shock,  with  the 
transformations 


r  =  R  +  v  s 
s 

(m-3) 

a  _  9 

Sr  v  Sb 


and  R  is  the  radius  of  the  shock  and  v  = 


oo 

"Ur" 


With  these  trans- 


s  p 

Moo  s 

formations  for  the  independent  variables,  we  can  now  consider  the 
derivatives  of  the  physical  quantities  v,  u,  etc.  in  the  stretched  radial 
direction  and  the  circumferential  direction  along  the  shock  have  the 
same  order  of  magnitude  as  explained  in  reference  13.  Substitute 
Eqs.  (IH-2)  and  (HI-3)  into  Eqs.  (DI-l),  equations  of  successive  orders 
can  be  obtained  consistently  with  respect  to  the  small  parameter  v.  In 
the  limit  v-*0>  we  obtain  the  differential  equations  for  the  leading  order 
shock  solution.  Expressed  in  the  nonstretched  coordinates  they  are  as 
follows: 


the  continuity  equation 


JZ  (pv)  =  0 


0H-4a) 


the  radial  direction  of  momentum  equation 


dp  ,  9v 
S7  +  pv  a~r 


0V  _  - 

—  W 


57  [(u'+  y  M)  |j]  » 


(m-4b) 
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the  circumferential  direction  of  momentum  equation 


8u 

0V  H7  - 


8 

Hr 


r 

[n-57] 


(m-4c) 


and  the  energy  equation 


9h. 

PV  87  = 


i £ 

8r 


^[k||]+(^+|,)[|X]a  +  ,[|£f  .  {m.4d) 


The  terms  in  Eqs.  (HI-4)  do  not  involve  the  derivatives  in  the  circum¬ 
ferential  direction,  however  the  dependent  variables  u,  v,  p,  p»T,  and  h 
are  functions  of  both  the  spacial  coordinates  r  and  0. 

Referring  to  Figure  1,  Eqs.  (IH-4)  can  be  integrated  along  a  ray 
of  constant  0  from  r  =  +  oo  to  a  point  r=  rx  (0),  to  be  defined  later. 

Eqs.  (IH-4a),  (IH-4b) ,  and  (H[-4c)  can  be  integrated  without  any  difficulty, 
since  both  the  left-hand  and  right-hand  sides  appear  in  the  forms  of  exact 
differential.  To  integrate  Eq.  (E!3-4d),  we  proceed  with  the  following 
rearrangement,  malting  use  of  Eqs.  (EH -4b)  and  (III-4c): 

Multiply  Eq.  (HI-4b)  by  v 


hence 


8p  ,  ,  t,  4  .  r8via  8  r/  . ,  4  .  8v,  8  .v3.  _  . 

V  8r  +  ^  +  3  ^  feP  =  87  ^  +  3^)v  87^  ~  pv  87  {T]  (m-5a) 


Multiply  Eq.  .(III-4c)  by  u 


pv  u 


8u 

H7  “ 


8u 


9ui 


r8u-i 


V  r  UUl  V  r  WUT  rvu-i: 

u  Hr  ^  Hr "  Hr  Hr^  " 


hence 


(m-5b) 


9h  9  pi  9T  *i,9  r  t  f  *  4  %  9v-i  ,  9  ■*  9u.^ 

pv  a?  =  a?  ^  W]  +  97  [{^  +  3  ^)v  a7]  +  87  [|JU  a7] 


8  rV2  n  o  ru~n 

"  Pv  37  L— ]  "  PV  [-7-] 


or 


8 

pv  87 


P's!  =  £  {*  ir - 


,  /,  4  .  8v  ,  8a  \ 
(U  +3  U)v  -57  +  ^ 


(m-4d,) 


u 


where  hg  is  the  stagnation  enthalpy  per  unit  mass  defined  as  hs  =  h  + 
Notice  that  daring  this  operation  no  postulation  has  been  made  on  the 
transport  coefficients  k,  y,  and  ju',  they  stand,  therefore,  as  from  the 
beginning,  as  the  LocaL  gas  properties  where  the  local  thermodynamic 
equilibrium  is  postulated.  Integrating  Eqs.  (HI-4a),  (IH-4b),  (HI-4c), 
and  (m-4d')  we  obtain  the  following: 


'+v3 


[pv]  =  p _ v  =  m 

LK  Jr=r1  ^00  00 


[p+mv-(n /+  |  m)  |^]r=ri  =  Pro  +  mv 


00 


8u 


mu_|J'87-‘r=r  =  mu 


00 


r  ,  ,  8T  ,  ,  4  .  8v  8ui  , 

[mhs"k  87  "  +  3  ^)v-§7  "  >itt,5FJ*=r1  =  ^ 


s  00 


(m-6a) 
(in- 6b) 

(m-6c) 

(m-6d) 


where  the  subscript  oo  in  the  right-hand  side  terms  denotes  that  the 
quantities  are  evaluated  at  r=  +  oo.  The  left-hand  side  terms  are  func¬ 
tions  of  r  and  0,  the  angle  of  inclination  of  the  ray.  Terms  such  as  v 
9v  ** 

and  -g-j-  now  can  be  viewed  as  independent  functions  of  0  if  the  integration 
limit  r.^  is  prescribed. 
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Inspecting  Eqs.  (HI- 6),  we  see  that  they  are  essentially  state¬ 
ments  of  balances  of  the  mass  flux,  the  components  of  stresses  and  the 
energy  flux  of  the  fluid  flowing  across  the  shock;  and  they  are  the  usual 
leading  order  shock  conditions.  Eqs.  (1H-6)  appear  essentially  the  same 
as  Eqs.  (5)  in  reference  14,  except  where  is  taken  as  zero  and  u=  u^ 
const,  across  the  shock,  which  is  the  condition  necessary  in  order  that  u 
shalL  be  bounded  as  we  integrate  along  the  stretched  radial  coordinate 
from  +  co  to  -co^^.  Therefore  should  be  retained  consistently  at 
r=  rx  while  ^  and  do  not  vanish.  Inspection  of  Eq.  (IH-6c)  shows 
that  the  implication  of  this  is  that  the  effect  of  change  in  the  tangential 
component  of  momentum  in  a  shock  does  not  come  in  the  same  order  as 
the  effect  of  the  shock  structure  (the  displacement  effect). 


It  has  been  pointed  out  in  the  Introduction  that  the  shock  detach¬ 
ment  distance  and  the  shock  curvature  are  different  parameters.  This 
means  that  the  shock  surface  near  the  axis  of  symmetry,  in  general, 
is  not  concentric  with  the  blunt  nosed  body  where  the  body  shape  near 
the  axis  can  also  be  approximated  by  a  sphere.  For  the  convenience  of 
study  the  flow  field  behind  the  shock,  system  of  Eqs.  (IH-6)  is  trans¬ 
formed  and  expressed  in  the  body-centered  spherical  polar  coordinate 
system,  where  a  parameter  characterizing  the  eccentricity  of  the  shock 
is  introduced.  This  is  explained  in  the  next  section. 


SECTION  IV 


COORDINATE  TRANSFORMATIONS 


Referring  to  Figure  1,  we  denote  the  distance  between  the  centers 
of  curvature  of  the  body  and  of  the  shock  on  the  axis  of  symmetry  by 
OgOfl=  d=  eRg,  where  Rg  is  the  radius  of  the  nose  on  the  axis  and  e  is 
the  dimensioniess  distance.  We  have  then  the  two  sphericaL  polar  coor¬ 
dinate  systems  (rfl,  0fl)  an-d  (rg,0g),  where  the  subscripts  indicate 
whether  the  system  is  shock-centered  or  body-centered,  respectively. 

From  the  geometry,  we  have  the  following  transformations: 


rg  =  [d2  +  rB3  +  2drBcos  0g] 


(IV -la) 


and 


sin  0s  =  rB  sin  0g[da  +  r03  +  2drBcos  9-g]"2 


rB  =  [d3  +  ra3-2drBcos  0j 


s  s 


sin  0  .  =  r  sin  0  [ds  +  r  3-2dr  cos  0  J 

ij  8  8  S  8  d 


(IV-lb) 

(IV  -2a) 
(IV -2b) 


The  radial  and  circumferential  components  of  velocity,  the  spacial 
derivatives  of  the  velocity  components  and  the  components  of  the  gradient 
of  a  scalar  quantity  at  any  point  P  can  be  transformed  from  one  system 
to  the  other  by  making  use  of  Eqs.  (IV-1)  and  (IV-2).  For  instance,  with 
the  velocity  components  defined  as 

d0g  _  d9B 

us  -  r8  ~ar  ub  rB  ~ar 

and 

drs  drB 

vs  -  "sr  vb  -  ~ar 
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the  circumferential,  component  of  velocity  referring  to  the  shock  system 
u8  is  then 

*  r.  3T 


,  8S.  “B  a8.  , 

=  r'[VB^i  ^  ^ 

vB[d  sin  0B]+uB[rB  +  d  cos  0B] 
[da  +  r_a  +  2dr„cos  0n]J 


da  dvg  dT 

Similarly,  the  quantities  vg  ,  -g-^  ,  -g-^—  ,  and  g-^  referring  to  the 

rs  s  s 

shock-centered  system  can  be  expressed  in  the  new  system  (rB»6g); 
they  are  Listed  in  Appendix  A.  The  parameter  d=  involved  in  the 

transformation  can  be  dtermined  uniquely  from  the  analysis,  and  is 
explained  in  the  later  part  of  the  work. 

The  system  of  shock  conservation  Eqs.  (m-6)  can  now  be 
expressed  in  the  new  coordinate  system,  with  terms  properly  defined 
as  referring  to  (rB»  Sg). 

The  continuity  equation  is 

p  v  =  -  F  = - e+(l+X)cos  6 - T  (IV -3a) 

[es  +  (l  +  y)s  +  2e(l+y)cos  0]2 

The  radial  component  of  momentum  equation  is 
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{— P  [hg-  -J  (va  +  u3)]  -  }  |e2  +  (l+y)2+2e(l+y)cos  0  j- 

=  [e+(l+y)cos  0]3  +  [e+(l  +  y)cos  0]  -^(1+y+e  cos  0)  v  -  (e  sin  0)  u  j- 

+  ea  -f(—-s—  [{e  sin  0)  v  +  (1  +y-e  cos  0)  u  ]  +  ~  [1+y+e  cos  0]2 

*  l  1+y  o y 

-  (e  sin  0)  [1+y  +  e  cos  0]+  (y^)  [easins8]-(^~-^)  ~  [1  +y+e  cos  0]  } 

(IV -3b) 


The  circumferential  component  of  momentum  equation  is 
[e+(l+y)cos  0]  ^[e  sin  0]v  +  [1+y+e  cos  0]  u  j- 

+  ei  {(±IT7^)  f(e  sin  “  “  (1+y+e  cos  0)  v] 

.  |v  _  (esinji)  [J+  cQ8  } 

1+y  90  1+y  1  1  J  90 

+  (e  sin  0)  [1+y+e  cos  0]  gy  +  [1+y+e  cos  0]2  -gy  j- 
=  (1+y)  sin  0  [e+(l+y)cos  0]  (IV-3c) 


The  energy  equation  is 


[e+(l+y)cos  0]  |es  +  (l+y)s  +  2e(l+y)cos  9  } {hs_ h8  qq  } 

e  r  9h  .  0  9h  -i 

+  (— — )  [ea  +  (l+y)3  +  2e(l+y)cos  0]  |(l+y+e  cos  0)  -(^-y^--)  ____  j 

rl  7  Y 
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+  (p— -)  {Pr  -1  j-  -£(l  +  y+e  cos  8)  v  -  (e  sin  0)  u  j- 


ra  '  3 


{* 


e  sin  Q)  j-^e  g)  v  +  (l+y_e  cos  8)  u  ]  +  [1+y+e  cos  0]s 

1  +y  ay 


-  ~(e  sin  8)  [1+y+e  cos  8]+  _  (g— 8-in  9)  [l  +  y+e  cos  0]  } 

8y  1+V  8  8  1+V  88  J 


+  y  80  1+y  80 


+  (p 1-)  -|p  -1  j-  -|(e  sin  6)  v  +  (1+y+e  cos  8)  u  j- 


ri 


j(®JLHLJ.)  [(e  sin  8)  u  -  (1+y+e  cos  0)  v]  -  (~ I) 

^  1+y  1+V  o0 


+y  30 


jejs^nJ5j  8^  [1+y+e  cos  gj  +  (e  g)  [i+y+e  cos  0] 


+  [1+y+e  cos  0]3  j-  =  0 


(IV -3d) 


ALL  the  bar  quantities  in  Eqs.  (IV-3)  refer  to  those  nondimen- 
sionaLized  with  respect  to  the  free  stream  quantities;  they  are  defined  as 


_  h  _ 

—  u  .  —  v  .  —  p  r~  s  7~ 
u  =  — -  ;  v=  —  ;  p  =  — k—  ;  h  =  —5  ;  h 


U 


U 


00 


and 


s  ua  s  00  u 


h  p _ 

s  00  —  00 

5—  >  P™  = 


r-RT 


y  = 


00  p  IT 
00 


e  = 


li  c 

T 


p  =  — 2  ;  P  =  1  P  ;  P  = 

r  A  4  lr 


LI  C 
^3  l 


«1  = 


;  ca 


Poo  PooURB 
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and 


M1  =  H  ;  Ha  =  U'  +  j  |i  - 

Again  Eqs.  (IV-3)  are  evaluated  at  y=y1(0)»  where  the  only  independent 
variable  is  0. 
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SECTION  V 


METHOD  OF  SOLUTION 


Using  the  boundary  layer  approximation  near  the  wall,  the  basic 
equations  of  single  species  Laminar  boundary  layer  flow  can  be 
expressed  as  folLows: 

Continuity  Equation 

-g—  (pv)  +  (p"u)  +  ~p  u  cot  0=0  (V-la) 

Momentum  Equation 

—  , —  8u  .  —  3u«  ^e  .  3  f—  3u.  ,,r  » 

P  (  u  ‘ST0  +  v  -0y)  -  “  -HT  +  e  ^  W  -By*  (V-lb) 


Energy  Equation 


p  (u 


^s 

ae 


+  V 


ah 
_ £ 

3y 


■)  =  e  | 


3_ 

l  3y 


(Ji- 

p 


3h 
_ e 

ay 


»  +  tv-1') 


3y 


Eqs.  (V-l)  are  the  same  as  those  appearing  in  reference  1  without  con¬ 
sidering  the  diffusion  phenomena,  except  here  all  the  dependent  and  the 
independent  variables  are  normalized  in  consistence  with  Eqs.  (TV-3). 

The  body  shape  is  approximated  by  a  sphere  with  radius  R„ ,  and 
r-RB 

y  =  -  .  The  problem  remains  now  to  seek  the  solution  of  Eqs.  (V-l) 

RB 

satisfying  the  usual  three  boundary  conditions,  the  temperature  distri¬ 
bution,  and  the  non-slip  boundary  conditions  at  the  wall,  and  the  four 

shock  conservation  Eqs.  (IV-3)  at  a  point  yet  to  be  determined.  Notice  « 

that  four  outer  boundary  conditions  appear  here  as  compared  with  the 

usual  three,  -  the  inviscid  veLocity  and  density  distribution  at  the  wall 

and  the  stagnation  enthaLpy^  prescribed  at  the  edge  of  a  thin  boundary 

layer.  This  extra  condition  suffices  to  determine  the  upper  integration 

limit  y  =  y1 ,  or  more  spectficaLLy,  the  detachment  distance  of  the  shock. 

The  solution  of  Eqs.  (V -1),  in  general,  can  be  obtained  by  machine 
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integration  with  suitabLe  transformation  for  the  dependent  and  the  inde¬ 
pendent  variables.  However,  the  problem  can  be  solved  more  easily  by 
using  a  technique  similar  to  that  employed  in  reference  4.  In  this 

technique,  the  physical  quantities  u,  v,  p,  h  and  their  radial  derivatives 

s 

can  be  expressed  in  terms  of  the  similarity  solutions  and  the  fictitious 
outer  flow  quantities  ug,  p  ,  and  hge  ,  which  in  turn  can  be  determined 
uniquely  from  the  present  analysis  (Figure  2). 

We  start  with  the  Lees'  transformation^ 


and 


puR^  sin  0  = 


9ty 

9y 


(V  -  2a) 


(V-2b) 


(V-2c) 


pvRg  sin  0  =  -  |-|  (V  -  2d) 

i  (0,y)  -J  2?  f(q)  (V-2e) 

hS/hse  =  .  <V-2f) 


where  quantities  with  subscript  e  refer  to  the  outer  flow  quantities  of 
the  boundary  layer  (i.  e.  ,  at  ti=  00)  with  the  boundary  conditions  as 
given  in  reference  1,  and  y  is  the  dimensionless  radial  coordinate. 
Consistent  with  Eqs.  (V-2),  any  physical  quantity  among  u,  v,  p  and  hfl 
and  their  spacial  derivatives  at  y  =  y1  {or  rj  =  T) x )  can  be  expressed  in 
terms  of  the  similarity  solutions  f  and  g  and  p  ,  u  ,  h  .  For  instance 

6  C  8  C 


r — I  —  df 

=  [  u  1  =  u 


y~y « 


dr\ 


TF  r\. 


=  u  f  * 
e  1 
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pflu,  rBu,  ueaPesin  9^ 

Vi-Vy=yi"  ^77fgi 


=  f/  [' 


9u 


-1  +  U  f  "■ 


Jk  {8in  0  pene} 


0  a  u  sin3 
Ke  e  e 


2? 


S 

-] 


re _ 

S  =  J  peueue  sina  0d0  » 

o 

where  the  bar  quantities  refer  to  the  quantities  nondimensionalized  with 
respect  to  the  free  stream  vaLues;  and  f x  '  and  f1"  denote,  respectively, 
first  and  second  derivatives  of  f  with  respect  to  r|  at  r|=  t)  .  Therefore, 
the  shock  conservation  Eqs.  (IV-3)  can  readily  be  transformed 
according  to  Eqs.  {V-2).  And,  together  with  Eq.  (V-2a),  this  forms  a 
complete  set  of  five  equations  for  the  new  variables  ue(0),  pe(9), 
hse(0)>  yx(0)>  and  1^(0). 

For  the  five  equations  to  hold  near  the  axis  of  symmetry,  the 
following  expansion  schemes  around  0=0  for  all  the  variables  can  be 
introduced  for  further  simplification  of  the  problem: 


ue(0)  = 

iue 

e  +  3u 

e  83 

(V-3a) 

Pe(0)  = 

ope 

+  ape 

e8 

(V-3b) 

5.e<®>  = 

e  +  a^se 

(V-3c) 

TljO)  = 

+  A 

e3 

(V-3d) 

7,(0)  = 

0*1 

+  A 

e3 

(V-3e) 

Introducing  Eqs.  (V-3)  into  the  transformed  shock  conservation  equations. 
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and  collecting  terms  of  equal  powers  of  0°  and  0a,  we  obtain  finally 
the  ten  algebraic  equations  for  the  ten  unknown  coefficients:  xue  ,  3ug  , 

0Pe  ’  3P  e'  <?se  *  a^se  '  <Ti  '  A  '  o*i  »  and  A  <or  e)*  They  are 


the  equations  on  the  axis  (coefficients  of  0°  terms) 

i 

[2e|j  p  u  l2  f  =  1 
L  oKe  l  ©J  l 


(V  -  4a) 


J  ^'gdri 


0^1  = 


2f. 


l  o^e  iue 


(V  -4b) 


(1+^)  0Pe]f/+  e1[1ue  opg]2  2f t  f"  ]  +  e  {[(xue  0Pe)f/- g] 


(V-4c) 


(  u  p  )f  7-g  f  g'  ^ 

+  {Pe  0J<J(2f1')(l+  -~)-(i^e  0fe>8(2W]  }  =  (l+<tfi>oP. 

1+0yi  Mi 


1  +  p  =  1=1  ^  h  +  Itl  [-fl]  +  e  u  [2f  '+2f  —) 

*00  r  oHe  o  se  2T  L  o  J  3  i  el  l  1  gx  (V-4d) 


r  o  se  ZT  1  a 

O ^ 


oKe 


_  _  _  _  2fxg/ 

[^s  oo  o^se  ®i  1  Ga  o^se  o^e  iue^  gt  ^  ’ 


(V -4e) 


and  the  equations  away  from  the  axis  (coefficients  of  03  terms) 


(V -5a-e) 


The  coefficients  K.  ,  ,  ©i  ,  R  ,  E.  (i=  1 , 2,  3, 4,  5)  of  the  linear 

Eqs.  (V-5)  are  functions  of  ^  ^  ,  0hBe  »  JV  ^  and  ^  (or  e): 

they  are  Listed  in  Appendix  B.  ALL  the  notations  appearing  here  are 
defined  in  consistence  with  Eqs.  (V -3) ,  (V^Z)^  and  (IV-3).  Notice  that 
in  Eqs.  (V -4)  and  (V -5),  the  quantity  e=  — ^  s  does  not  appear  as  an 
extra  variabLe  other  than  .  Since,  from  t^e  matching  surface 

_  *,(0) 

ri(e)  *  -V-  = 1  +  y.<®>  =(1  Vi  > +  *  ®2 

B 

and 

(1+  y  +e)  =  (shock  radius  at  0=0)  = 
o  i 


it  can  easily  be  derived  that 


_  dr  _  d3r 


1 0=0  , 
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Inspection  of  Eqs.  (V-4)  and  (V-5)  shows  that  the  two  sets  of 

equations  are  coupled  in  a  sense  that  Eqs.  (V-4)  involve  the  six 

variables  „u  ,  p  ,  h  ,  ri  ,  _y  and  e  while  Eqs.  (V-5)  involve  all 

l  e  oHe  '  o  se  o  h  o7i  _  ^ 

■  o  g  2  P  g  2  h  q  p 

the  ten  variables  in  which  — rr  -  ,  — ,  — == —  ,  and  ri.  appear  linearly. 

l  ue  oPe  ohse  2  1 

This  suggests  the  following  iteration  procedure  for  the  solution  of  the 
problem:  for  prescribed  initial  value  of  e,  values  of  u  ,  p  .  h  , 

Qr)1  ,  and  Qy  can  be  determined  by  solving  Eqs.  (V-4);  then  the  new  value 
of  e  can  be  obtained  by  solving  one  of  the  five  equations  in  Eqs.  (V-5)  for 
the  next  cycle  of  iteration.  The  numerical  example  i«  presented  in  the 
next  section. 
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Inspection  of  Eqs.  (V-4)  and  (V-5)  shows  that  the  two  sets  of 

equations  are  coupled  in  a  sense  that  Eqs.  (V-4)  involve  the  six 

variables  u  ,  p  ,  h  ,  ri  ,  y  and  e  while  Eqs.  (V-5)  involve  all 
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the  ten  variables  in  which  3_—  ,  --P_e-  ,  ,  and  rj  appear  linearly. 

iue  o  Pe  o^se 

This  suggests  the  following  iteration  procedure  for  the  solution  of  the 
problem:  for  prescribed  initial  value  of  e,  values  of  ug  ,  0pg  ,  Qhse  » 

Qr)1  ,  and  can  be  determined  by  solving  Eqs.  (V-4);  then  the  new  value 
of  e  can  be  obtained  by  solving  one  of  the  five  equations  in  Eqs.  (V-5)  for 
the  next  cycle  of  iteration.  The  numerical  example  i:s  presented  in  the 
next  section. 
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SECTION  VI 


NUMERICAL  EXAMPLE  AND  DISCUSSION  OF  RESULTS 

Eqs.  (V-4)  and  (V-5)  are  applied  in  the  specific  calculation  for  the 
case  of  Mqo=  5.7  and  Tq=  free  stream  stagnation  temperature  =  2100°R. 

In  performing  the  computation,  the  cold  wall  approximation  has  been  used 
for  the  enthalpy  profile  g  =  f and  numerical  values  for  the  Blasius 
velocity  profile  f  have  been  obtained  from  reference  17.  The  linear 
temperature -coefficient  of  viscosity  law  has  been  assumed.  {Notice  that 
the  use  of  Blasius  velocity  profile  does  not  imply  the  present  scheme  is 
restricted  only  to  zero  tangential  pressure  gradient.  In  fact  the  present 
analysis  is  completely  general  for  various  tangential  pressure  gradient 
imbedded  in  the  momentum  equation  of  the  boundary  layer  equations.) 

The  results  are  presented  in  Tables  I  and  II  for  the  range  of 
Reynolds  number  Res  =  .Rs'^^s  oo  (reference  4)  from  20  to  400. 

_  Ht 

Where  the  and  the  Qh  are,  respectively,  the  density  and  the  stag¬ 
nation  enthalpy  at  the  matching  point  on  the  axis.  They  are  plotted  as 
functions  of  the  Reynolds  number  in  Figure  3.  It  is  interesting  to 
observe  that  the  value  of  p  decreases  as  Re„  increases  and  that  it 
approaches  the  value  corresponding  to  that  across  a  discontinuous  shock 
(^  =  6.22626  based  on  Moo  =  5.7  and  F=  1.325  after  the  shock).  The  value 
of  0p  increases  rapidly  as  the  Reynolds  number  decreases;  the  unbounded 
trend  is  due  to  the  fact  that  the  cold  wall  approximation  has  been  used  in 
the  calculation  and  ^  0pe/f The  stagnation  enthalpy  at  the 

matching  point  0hsl  increases  as  the  Reynolds  number  increases,  and 
approaches  the  value  of  the  free  stream  stagnation  enthalpy.  This  indi¬ 
cates  the  fact  that  the  heat  flux  in  the  radial  direction  after  the  shock 
becomes  more  and  more  important  as  the  Reynolds  number  decreases. 

The  results  of  p  /  p  and  h  /  h  are  also  shown  in  Table  I.  The 
negative  value  indicates  both  the  stagnation  enthalpy  and  density  decay 
along  the  circumferential  direction.  The  rapidness  of  the  decay  can 
also  be  observed  from  the  magnitude  of  the  coefficients.  As  both 


3fh/ofh  an<^  s^se/o^se  approach  zero  as  the  ReynoLds  number  increases, 
one  can  easily  observe  that  the  concept  of  the  constant  inviscid  density 
distribution  is  reasonably  justified  at  high  Mach  number  and  high 
Reynolds  number.  The  shock  eccentricity  e  is  plotted  against  r|j 
as  shown  in  Figure  4.  The  matching  distance  Qy1  on  the  axis  is  plotted 
agains  Res  as  shown  in  Figure  5.  The  value  of  Qy  stays  within  3  per¬ 
cent  for  the  range  of  Res  from  50  to  400,  and  is  of  the  order  of  one- 
half  of  the  inviscid  detachment  distance.  This  can  probably  be  explained 
as  the  shock  is  merging  into  the  boundary  layer.  The  accuracy  of  this 
result  has  yet  to  be  verified  with  the  experimental  data  at  this  Mach 
number.  The  rapid  decrease  of  at  Reg  =  20  indicates  the  probable 
range  of  application  of  the  present  scheme,  as  it  seems  to  be  unrealistic 
to  extend  the  range  of  continuum  approach  any  further. 

In  terms  of  the  boundary  layer  language,  the  surface  heat  transfer 
value  is^ 


nm  k  k  9h 

a  =  k  JlL}  -  w  oh.  _  w  _ s. 

W  3yW  c  8yW  "  c  9y  W 

Pw  Pw 


k  „  knur 

w  ,  E  w  Hw  e  o 


h.eS'<°>Ow=  - 


w 


w 


- I  hse8/(0) 

(2s)2  se 


For  the  same  wall  temperature,  the  heat  transfer  ratio  based  on  the 
boundary  layer  value  at  the  stagnation  point  can  easily  be  derived.  It  is 

»  o  u  i  h 

_  =  &J2)  (L-Sj2  £■ *?.) 

^B.L.  o^s  iUo  ^s  oo 

where  Qp  '  is  the  normalized  inviscid  stagnation  density  after  the  dis¬ 
continuous  shock,  and  ^I0  is  the  inviscid  velocity  distribution  based  on 

reference  1.  For  M  =5.7  and  P=  1.325,  we  have 

oo 


p'  =  6.67518  u  =  0.52123  . 

Ob  1  O 


The  resuLt  of  -■ — 3 —  j.s  shown  in  Table  n.  Figure  4.1  of  reference  15 
9B.L. 

is  reproduced  in  Figure  6  for  comparison  of  the  present  theory  with 
the  theory  developed  in  reference  12,  as  well  as  with  the  experimental 
data  of  references  4,  18,  and  19.  The  plot  shows  a  slight  improvement 
on  the  heat  transfer  than  the  theory  developed  in  reference  12  in  the 
range  of  application  of  the  present  scheme.  This  indicates  that  the 
shock  eccentricity  does  not  play  a  very  important  role  as  far  as  heat 
transfer  value  is  concerned.  The  theory  breaks  down  beyond  the  point 
where  Res  is  of  the  order  of  400,  as  apparently  one  cannot  expect  the 
"direct  matching"  to  be  a  realistic  scheme. 


The  heat  transfer  value  is  also  converted  and  expressed  in 
Cha/Rcs  ,  where  the  dimensionless  heat  transfer  value  Cjj  = 

The  result  is  plotted  against  Res  as  shown  in  Figure  7.  Some  of  the 
experimental  results  of  references  4,  18,  19,  and  20  are  also  included 
for  comparison.  It  is  seen  that  at  the  large  Reynolds  number  end,  the 
difference  on  the  heat  transfer  value  between  the  theories  of  references 
4  and  14  is  only  of  the  order  of  two  percent  when  the  absolute  scale  is 
used. 
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SECTION  vn 


CONCLUSION 

It  was  found  that  the  heat  fLuxes  and  stress  components  behind 
the  shock  have  strong  influences  on  the  flow  field  around  the  body  in 
the  range  of  Reynolds  number  of  the  present  study.  The  heat  transfer 
value  decreases  monotonically  towards  the  free  molecular  value  as  the 
Reynolds  number  decreases.  The  shock  eccentricity  was  found  to  have 
very  little  influence  on  the  heat  transfer  value.  The  flow  density  at  the 
matching  surface  increases  rapidly  as  the  Reynolds  number  decreases. 
On  the  other  hand,  the  stagnation  enthalpy  at  the  matching  surface 
decreases  as  the  Reynolds  number  decreases.  Both  the  stagnation 
enthalpy  and  the  density  at  the  matching  surface  yield  the  regular 
Rankine-Hugoniot  values  at  the  large  Reynolds  number  end. 
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APPENDIX  A 


TRANSFORMATIONS  BETWEEN 
SHOCK-CENTERED  SPHERICAL  COORDINATE  SYSTEM 
AND  BODY-CENTERED  SPHERICAL  COORDINATE  SYSTEM 


1.  TRANSFORMATIONS  OF  THE  SPACIAL  COORDINATES 

From  the  geometry  (referring  to  Figure  1),  the  coordinate 
transformations  between  the  two  systems  (rg,  0  )  and(rB>0B)  are 

(n-la) 

(n-ib) 

(II -2a) 

(H-  2b) 


r  =  [d®  +  r  3  +  2dr  CCS  0-n  ]2 

S  ±J  ±5  X5 

sin  6g  =  rB  sin  [d®  +  +  2drB  cos  Bg]”2 

and 

l 

rn  =  [d2  +  r  3  -  2dr  cos  0  ]2 
B  1  s  s  s 1 

sin  0„  =  r  sin  0  [da+r  3-2dr  cos  0  1  2 
B  s  sL  s  s  sJ 


where  d  is  the  distance  between  the  centers  0._  and  O  ,  and  the  subscripts 

.D  S 

B  and  s  refer  to  the  coordinate  system  and  are  body-centered 
or  shock-centered,  respectively. 


2.  TRANSFORMAT  IONS  OF  THE  VELOCITY  COMPONENTS 


The  radial  and  circumferential  velocity  components  are  defined  as 


u 


s 


r 

s 


Ut 


rB 


d6B 

dt 


31 


and 


Then, 


dr 

s 

v  =  - 

s  dt 


= 


drB 

dt 


Us  rs  dt  ^MrB'6B^ 


90„  drT 


r.  fc 


9  0„  deT 


B 


dt 


aeT 


dt 


9e  90 

r  8  ,  d  S  i 


vfi[d  sin  0  ]  +  u  [r  +d  cos  0  ] 

Ug  = - y - 

[da+rBa  +  2drB  cos  0B  ]2  (A-l) 

Similarly,  we  obtain 

vR[rR+d  cos  0_]  -  u-r,  [d  sin  0_] 
vs  =  -5-5 - 5 - n - 1_5_  (A-  2) 

[ds+rB2+2drB  cos  @B]2 

3.  TRANSFORMATIONS  OF  THE  DERIVATIVES  OF  THE  VELOCITY 
COMPONENTS 


9u  dQ  9u  9r_  9u 

__i  =  (_B )  _JL  +  (_B )  _J 

9r  '9r  90^  V9r  8rT 

s  s  B  si 
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> 


-  d  sin  0 


-  d  sin  0T 


90r 

■P  _ _ 0 _  _  _  JL J 

8rs  [rs3  +  d3-2drs  cos  0g]  rB[rBa  +  ds+ 2drB  cos  0B]2 


and 


8r 


B 


rB+d  cos  0B 


9rs  [r_s  +  da  +  2drT,  cos  0„]2 


Making  use  of  Eq.  (A-l),  we  obtain 


9u 


d  sin  0, 


9r  r„ Jr^+d3*-  2drB  cos  0Q) 


[(d  sin  6B)uB~(rB+d  cos  SgWg] 


s  B'  B 


3u. 


+  [rB3  +  d3+2drB  cos  0B]-1  {|>B  +  d  cos  @B]3 


dvT 


+  (d  sin  0Q)  [rB  +  d  cos  0B] 


d  sin  0  8 Up.  d3sin30  3v  . 

-  ( - [r  +d  cos  0  ]  —2.  .  ( - B)  B  | 

r-o  B  B  80  r  8@  J 


Similarly,  we  obtain 


(A- 3) 


8v 


d  sin  0T 


9rs  rB(rBs  +  ds+ 2drBcos  @B) 


[(d  sin  @B)vB+(rB-a  cos  Ogjug] 


3vt 


+  [rB3  +  d3+2drBcos  0B]_1  {[rfi+d  cos  @B]3 

B 
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da  sin3  0_  9uT 


+  (- 


-) 


lB 


9  6t 


d  sin  0_  9v_  . 

'  < - 7—)  [rB  +  d  cos  aeT  } 

B  B 


(A-4) 


4.  TRANSFORMATION  OF  THE  DERIVATIVE  OF  TEMPERATURE  IN 
THE  RADIAL  DIRECTION 


Thus, 


9T 

9r 


86U 

=  fc-^) 


9T 


9r  9  0, 


8rT 


9T 


9r  9r, 


9T 


[rB+d  cos  0B] 


9T 


9rs  [ 3 -f da  +  2drB  cos  0g]2  9rB 


[d  sin  0g] 

[rB3  +  d3+2drB  cos  0B P  rg 


9T 

"b 


(A- 5) 


With  Eqs.  (A-l)-(A-5),  the  shock  conservation  Eqs.  (1-6)  in 
Section  I  can  readiLy  be  expressed  in  the  (rB>  @B)  coordinate  system 
as  in  the  form  of  Eqs.  (II-3). 
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APPENDIX  B 


COEFFICIENTS  OF  THE 
SHOCK  CONSERVATION  EQUATIONS  (EQS.  Y-5) 


In  the  following  expressions  for  the  coefficients,  aLl  the  terms 
invoLve  f  and  g  and  their  derivatives  are  evaluated  at 
ri=  0iv  and  3y1=  f 

o  7  1 


K^=^[2  o^x^-13 


Ka  =  3^-  [2oT1l(^-)-l] 
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TABLE  I 


DENSITY  AND  STAGNATION  ENTHALPY  AS  FUNCTIONS 
OF  THE  MATCHING  DISTANCE  (BOUNDARY  LAYER  COORDINATE) 
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FIG.  2.  Diagrammatical  Drawing  of  the  Matching  Velocity  Profiles 
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FIG.  3.  Density  and  Stagnation  Enthalpy  at  the  Matching  Point 
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Res=  REYNOLDS  NUMBER  BASED  ON  REF.  4. 
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100  1000 
FIG.  7.  Heat  Transfer  Value  Cjj/Re^"  Versus  Res 
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